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Abstract 

We study the s-channel production of a single top squark in hadron colli- 
sions through an i2-parity-violating mechanism, examining in detail the case 
in which the squark decays through an imparity-conserving process into a 
bottom quark, a lepton, and missing energy. We show that the top squark 
can be discovered if its mass is less than 400 GeV, or that the current bound 
on the size of the i?-parity-violating couplings can be reduced by up to one 
order of magnitude with existing data and by two orders of magnitude at the 
forthcoming run II of the Fermilab Tevatron. 

PACS numbers: 11.30.Pb, 12.60.Jv, 14.80.Ly 
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I. INTRODUCTION 



In supersymmetric extensions of the standard model, particles may be assigned a new 
quantum number called R parity (R p ) (]J . The particles of the standard model are Rp even, 
and their corresponding superpartners are R p odd. If R p is conserved, as is often assumed, 
superpartners must be produced in pairs, each of which decays to a final state that includes 
at least one stable lightest supersymmetric particle (LSP). If R p is not conserved, then 
many of the experimental signatures of supersymmetry that are usually studied may not 
be observable. Hence, it is important to determine whether R p violation occurs and, if so, 
whether the effects are significant in collider experiments. Current bounds on possible R p - 
violating couplings, obtained principally from quantum corrections to particle decays and 
neutral meson mixing, are relatively restrictive for the first generation of quarks and leptons 
but are much less so for states of the second and third generations @||. 

Searches for R p violation often focus on the decays of superpartners under the assumption 
that they are produced in pairs via an i? p -conserving interaction 0. However, squarks and 
gluinos are generally predicted to be relatively heavy, and therefore their pair-production 
cross sections incur large phase space suppressions at Fermilab Tevatron energies. 

In this paper, we expand and improve upon our previous proposal Q to probe R p - 
violating couplings by searching for the s-channel production of a single squark through 
an R p - violating mechanism ||. The greater phase space afforded by the production of one 
heavy particle, instead of two, will allow searches to reach much higher masses than those 
based on traditional squark pair-production. Further, since the i?-parity violation occurs 
in production, the cross section will provide a direct measure of the individual _R p -violating 
couplings, rather than of products of different _R p -violating couplings. 

In Sec. |Tj, we describe the baryon-number-violating portion of the i? p -violating La- 
grangian and establish our notation and conventions. The cross section for pair-production 
of top squarks is small compared to that for the process qq' — > t\ that can proceed through 
.Rp-violating couplings. In order to capitalize on the enhanced rate for single production 
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it is necessary to examine the decay products of the top squarks. We address in Sec. |T| 
the possibility that decay will occur through an ^-violating process into two hadronic jets 
and show that this rate is far below that for standard quantum chromodynamics (QCD) 
production of jet pairs. 

In Sec. [TV], we study observability of single top squarks by focusing on the clean R p - 
conserving decay, t\ — > bxf, with xf — > l + + v i + Xv Here, / is an electron or muon, and the 
Xi and Xi are the lowest-mass chargino and neutralino states of the minimal supersymmet- 
ric standard model (MSSM) ||. We introduce a method to extract a top-squark mass peak 
from the standard model background, principally W+ hadronic jets. To study the effects of 
parton showering, hadronization, and color-recombination, we simulate the signal with the 
HERWIG [0] Monte Carlo event generation program. The response of a typical experimen- 
tal detector is modeled using the full SHW detector simulation package ||. We explore a 
wide range of the parameter space of the minimal supergravity (MSUGRA) model of super- 
symmetry breaking |J. Each of these improvements renders the analysis more realistic and 
strengthens the conclusions reached in our earlier work [|J. 

Our conclusions are stated in Sec. |V|. For top squarks with mass mr < 400 GeV, we 
show that discovery at the level of five standard deviations (5a) is possible at v^S* = 2 TeV 
with an integrated luminosity of 2 fb _1 , provided that the R p - violating couplings A" > 0.02 
- 0.1. Otherwise, a 95% confidence-level exclusion limit can be set for A" as a function of the 
mass up to about 500 GeV. For the lower integrated luminosity and energy of the existing 
run I data, values of A" > 0.05-0.2 can be excluded at the 95% confidence level if m? = 
165-350 GeV. These limits would constitute a significant improvement over the current 
95% confidence-level upper bounds of A" < 1 for squarks of the third generation 
With the large increase in cross section at the CERN Large Hadron Collider (LHC), the 
discovery region and limits may be improved by another order of magnitude in A" and to 
top-squark masses approaching 1 TeV. The discovery reaches, or exclusion limits, at the 
Tevatron and LHC depend on the value of the top-squark mass, largely independent of 
other supersymmetric parameters. 
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A lower limit of about 165 GeV on the accessible value of the top-squark mass via our 
method is determined by the selections made on the b jet and lepton momenta in order to 
achieve a satisfactory signal to background ratio. An upper bound of about 550 GeV is set 
by the event rate. To explore values of the top-squark mass well below that of the top quark, 
it appears necessary to combine results of our analysis with those based on pair production 
of top squarks. 

For the values of the i?-parity-violating couplings of interest to us, the i?-parity- violating 
width of the Xi remains negligible and for all practical purposes the Xi LSP is stable. We 
leave to future work the exploration of larger values of i?-parity-violating couplings and 
consequent decays of the x± within the fiducial region of detectors. 



II. MSSM SUPERPOTENTIAL AND NOTATION 

In general it is possible to have i? p -violating contributions to the MSSM superpotential 
that violate baryon- or lepton-number. However, limits on the proton decay rate severely 
restrict their simultaneous presence. In this paper, we assume the existence of only a baryon- 



number-violating coupling of the form |10 



Wfe = KkU'Dpt . (1) 

Here, Uf and D\ are right-handed-quark singlet chiral superfields, i,j, k are generation 
indices, and c denotes charge conjugation. This form of the superpotential fixes our choice 
of normalization of X"j k . It is equal to where A^- is the coupling used in HERWIG 0. 

In four-component Dirac notation, the Lagrangian that follows from this superpotential 
term is 



+ d Rkl v^P R d jP + /i.e. , (2) 



where j < k. 
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For production of a right-handed top squark via an s-channel subprocess djdk — » 
the relevant couplings are X r ^ 12 , Ag 13 , and A 323 . The most direct limits on these couplings 
come from the measurement of Ri, the partial decay width of the Z boson to hadrons over 
its partial decay width to leptons. The analysis of Ri provides 95% confidence-level upper 
bounds of \'^ k < 1 HH. The Bayesian limits^] given as a function of mass in Ref. |§, 
are approximately A 3 ' 12 < 1.52 + 0.425 x (m^/100 GeV), and A 3 ' 13 323 < 0.585 + 0.162 x 
(rr^/100 GeV). 

Subsequent to our Letter |4j], it was argued that rare B + decays and K°-~K° mixing 
constrain A 3 fe strongly ||11|| . The analysis of these decays indirectly constrains products of 



different A^, whereas the method we propose may be used to directly probe each coupling 
independently. 

Shown in Fig. [I] is the Feynman diagram for s-channel production of a single top squark 
via the partonic subprocess d + s — > £r. The right-handed top-squark interaction state is 
related to the mass eigenstates by tn = —ti sin 9j + t 2 cos 9j. For simplicity of notation, and 
motivated by models where sin% pa 1, we present results in terms of a measurement of t\. 
In general, the results are valid for whichever mass eigenstates contain some amount of the 
right-handed top squark. 

The color- and spin-averaged cross section for inclusive t\ production is 

^ = §sin 2 %^|A 3 ^|% fc (r), (3) 

j<k 

where S is the square of the hadronic center of mass energy, and r = m~ / S. The integrated 
parton luminosity $>jk contains a convolution of the parton distribution functions (PDF's): 
g? Cg) s, d®b, and s <E> b, where d, s, and b label the PDF's of the down, strange, and bottom 



Ref. g, it is stated that a classical statistical analysis rules out the X'^ 3 couplings at the 2a 
level, and places a bound of A 3 ' 12 < 2.7. However, we do not agree that the analysis of Ref. [fj 
places a significant bound on A^-g. 
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quarks, respectively. The cross section for t 2 production is the same as above but with sm 2 6j 
replaced by cos 2 % 

The mass dependences of the next-to-leading order (NLO) cross sections for i? p -violating 
single production ]n| and i? p -conserving pair production [TJJ of top squarks differ signifi- 
cantly, as shown in Fig. |]. Results are presented for the sum of t\ and t\ production and 
couplings \'^ k = 0.1. Even if X'^ k is reduced to 0.01, two orders of magnitude below the 
current bound, the single-top-squark rate exceeds the pair rate for mj > 100 GeV. The 
parton luminosities determine that the contribution to the total cross section of the terms 
proportional to Ag' 12 : Ag 13 : A£ 23 is about 0.72 : 0.24 : 0.04 at m~ ti = 200 GeV. 

In order to simplify the explanation of some of the general features of this study, we 
use the notation A" = A 3 ' 12 = X^i3 = ^323- However, our method is identical even if only 
one of the three couplings is non-zero. The results for independent couplings are related 
by a simple functional form and are presented below. Since the cross section in Eq. ([$]) is 
proportional to the square of the product \'^ k sin % a model-independent analysis of the 
data limits the product. Correspondingly, limits for £ 2 are really on \'^ k cos 9j. In Sec. fVJ we 
present minimal supergravity models of supersymmetry breaking in which the mixing angle 
9j is uniquely determined by the model parameters. Within the context of these models, 
limits are placed on the actual X'^ k couplings. 

III. i?-PARITY- VIOLATING DECAYS 

In the .Rp-violating MSSM, the right-handed up-type squark ur can decay into quark 
pairs ur — > dj + dk via the A" couplings. The decay width of the top-squark mass eigenstate 
into two jets is 

r(tx-ii) = gsin 2 ^|Ay 2 . (4) 

j<k 

If A" > 2, the width of the top squark is greater that its mass. This is true generally for any 
squark and any baryon-number-violating coupling. Hence, there is a practical upper limit 
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on A" for the model to contain narrow resonances. 

In Fig. ^, we present the dijet mass distribution that results from the i? p -violating decays 
of singly-produced top squarks with A" = 1. We simulate the signal for various masses with 



a matrix-element calculation that uses HELAS |14| subroutines. For comparison, we also 
reproduce data published by the Collider Detector at Fermilab (CDF) Collaboration |T5 . 
We apply the same cuts in our simulation as were used on the data. Figure [| shows that the 
signal to background ratio (S/B) in the dijet channel is less than 1/100 for any mass. For 
A" different from 1, this ratio is reduced further to S/B < 0.01 x (A") 4 / [(A") 2 + c], where 
c is a constant proportional to the ^-conserving width. Because of the overwhelming jet 
rate from standard strong interaction processes, we confirm the expectation || that the 
dijet data can neither exclude i?-parity-violating production and decay of top squarks nor 
provide useful bounds on A". 

IV. i?-PARITY-CONSERVING DECAYS 

In this section we focus on probes of A" couplings through i? p -conserving decay modes 
of the top squark. If kinematically allowed, a heavy top squark decays into a bottom quark 
(b) and the lightest chargino xt ■ If m t 1 < m x + + m &> the decay would occur into a charm 
quark and a neutralino via flavor changing loop diagrams []nj. For the top-squark masses 
of interest to us, and in the context of the MSUGRA model in which we work, the decay 
t\ — > b + xt dominates. Other channels open at large t% masses (e.g. decay into a top 
quark and the lightest neutralino, t + x\)i but their decay widths are less than half of the 
R p - conserving decay width for the masses we consider. The chargino typically undergoes 
a three-body decay into a neutralino and the decay products of a W boson. Since it is 
challenging to observe an all-jets plus missing energy mode, we concentrate on the case 
where the final state is a bottom quark plus a neutralino and an electron or muon from the 
chargino decay. 

We show the branching fractions of the top squark into two jets and of the cascade 
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(ti — > bxl — » bl + viXi) as a function of A" in Fig. |j. The i? p -violating branching fraction of 
the chargino is also shown. The i? p -conserving decay dominates for small A". For both the 
top squark and the chargino, the width into jets increases when A" is large. The branching 
fraction into the signal mode is inversely proportional to (A") 2 . However, as seen in Eq. 
the production cross section is proportional to (A") 2 , and thus the entire cross section behaves 
as 

a{pp -> bl + vo$) = a(pp -> t x ) x BR(ti -> &xl) 

~ A A" 2 + ar Kp A" 2 + &rf* ' ^ 

In Eq. (D, Y~ p and are the R p - conserving partial decay widths of the lightest top squark 
and chargino, respectively, and a and b are appropriate proportionality factors. 

In the presence of _R p -violating couplings, the neutralino is not stable. There are three 
decay scenarios that potentially give rise to distinct signatures in the detector. In the first 
scenario, only the A3 - fc couplings are larger than 10~ 4 -10 -5 . In this case, the neutralino must 
decay through an off-shell top quark and off-shell top squark into a five-body final state. A 
neutralino whose mass is 100 GeV or less has decay lifetime cr > 100/(A") 2 m (see Fig. 26 
of the first paper of Ref. 0) and decays far outside the detector. Hence, the signature of 
the neutralino is simply missing energy in the detector. In this paper we address the case 
where the neutralino either decays outside of the detector volume or leaves remnants that 
are too soft to identify. The signal is a tagged 6-jet, an electron or muon, and large missing 
transverse energy. The main backgrounds are single-top- quark production, and H^+jets 
processes that either contain a bottom quark or have a jet that is mis-tagged as a 6-jet. 

If other baryon- number- violating couplings are large, then the neutralino can decay inside 
the detector. The signature of this second scenario includes two or three extra jets that 
appear to come from the primary vertex. The modeling of the energy distributions of these 
jets depends on explicit choices of some parameters. In an interesting third possibility, at 
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least one of the first- or second-generation couplings, A" 23 , A^, is greater than 10~ 3 . The 
decay may either have a displaced vertex or occur in a calorimeter. These possibilities 
warrant further study. 



A. Simulation 

In order to include the effects of parton showering and hadronization, we simulate the 
signal with the Monte Carlo program HERWIG 6.1 0, a version that includes _R p -violating 
interactions, and we use the SHW 2.3 |§ detector simulation package. The SHW package 
provides a reasonably close approximation of the expected acceptance of the upgraded CDF 
and DO detectors to signal and background processes. It determines what charged tracks and 
calorimeter energies the detector would record and supplies information about the trigger 
and reconstructed states such as electrons, muons, and hadronic jets, including b- and c-jet 
tagging. 

In addition to interfacing SHW to HERWIG, we make a few modifications for this 



study. We include the NLO ^-factors calculated in Ref. [|T^]. A kx cluster algorithm [17 
for hadronic jets is added to provide an infrared-safe way of combining calorimeter towers. 
We use a k? cone size of 1, similar to a fixed cone size of AR = 0.7. We also correct for 
the typical out-of-cone and threshold energy losses in the jet reconstruction of about 4 GeV 
per jet. After this correction, we find that the radiation, showering, and color-reconnection 
modeled in HERWIG 6.1 have little effect on the distributions of the reconstructed objects 
with respect to the parton level results presented in Ref. [|J. 

We use impact-parameter tagging as defined in SHW to tag 6-jets and to acquire the 
energy-dependent rate for charm quarks to fake bottom quarks. For a conservative estimate 
of the background, we choose the mis-tag rate for light-quark jets to be the greater of 0.5% or 
the output of SHW. The event trigger is either an electron or muon that passes the SHW 
triggering condition listed in Table |, typically a transverse energy Eti > 15 GeV. The 
missing transverse energy E T is calculated for the entire detector (calorimeters plus muon 
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chambers) after energy corrections. In order to stay above the level of detector fluctuations 
we choose E T > 20 GeV. The final significance is insensitive to the jB T cut as long as it is 
not raised above about 0.4 x mi. 

In Table |, we provide the geometric acceptance of a detector as defined by SHW. We 
also list the minimum transverse energy used to define identifiable objects (6-jets, leptons, 
or "hard" jets). The lepton-finding algorithm in SHW is based on the CDF run I geometry 
and thus covers a slightly smaller pseudorapidity region than used in our original Letter 
(see Table I of Ref. [|J). There is an additional loss of 20% of the leptons compared with 
our initial assumptions due to modeling of detector efficiencies. After these detector effects 
and the slightly different choice of cuts are accounted for, the full simulation is in complete 
agreement with the exact leading-order matrix element calculation used in our Letter [[|. 
We also confirm that adding Eq. (|3|) to PYTHIA 6.131 |Tj| produces the same results, an 
independent confirmation that color-recombination ambiguities, present for finite baryon- 



number- violating couplings ||19|| , do not impact the results. 

The backgrounds are modeled with tree-level matrix elements obtained from MAD- 
GRAPH |2Uj convolved with leading-order CTEQ5L |2T|] parton distribution functions, at a 



hard scattering scale fi 2 = s. In order of importance, these backgrounds arise from produc- 
tion and decay of the standard model processes Wc, with a charm quark c that is mistaken 
for a b; Wj, with a hadronic jet that mimics a b; Wbb; Wcc; and single-top- quark production 
via Wg fusion. The tagging efficiencies and geometry are the same as used in SHW. In an 
experimental analysis, the Wj background would be normalized by the data. To simulate 
the resolution of the hadron calorimeter for background events, we smear the jet energies 
with a Gaussian whose width is AEj/Ej = 0.80/ ^^Ej © 0.05 (added in quadrature). To 
verify the accuracy of our background estimation, we perform a full HERWIG and SHW 
simulation for the Wj background. The relevant distributions match those of the matrix el- 
ement calculation to better than 10%. As explained below, numerical limits on the couplings 
are insensitive to small uncertainties in the background. 

The b quark recoils against the chargino in the primary decay of the top squark. Thus 
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the measured Et spectrum for the 6-jet is peaked near the kinematic limit 

?tt| — m?. + + ml 

7-imax _ *l Xi //?\ 

Et » ~ 2^ ■ (6) 

An estimate of the mass difference between the top squark and chargino may be obtained 
if a prominent peak is found in the Exb spectrum. If we invert this equation, a kinematic 
constraint appears on the lowest top-squark mass that may be probed with this method, 
m~ tx > m~+ + Efp£. 

In Fig. ||, we show the E T spectrum for the 6-jet, tagged by the SHW impact-parameter 
method, for two different top squark masses, as well as the Wbb background. The background 
falls rapidly as a function of Exb- The 6-jet from the single-top-quark background peaks near 
60 GeV, but it makes a small contribution after cuts. We choose a cut of 40 GeV on the E? 
of the 6-jet to ensure a reasonable tagging efficiency (greater than 50%). Variation of this 
number between 30 GeV and 80 GeV has no effect on the significance of the signal as long 
as the signal remains kinematically allowed. 

The lepton from chargino decay (xt ~^ ^ +1J iXi) tends to be soft since much of the energy 
goes into the mass of the neutralino, and the rest is split by the three-body decay. In 
contrast, the lepton in the background comes from decay of a real W. It clusters around 40 
GeV (1/2 of the W mass) and has a broad rapidity distribution. After reconstruction, the 
remaining events are distributed from 20 GeV to 70 GeV (see Fig. |6]). Given the different 
shapes of the pr spectra, we impose a hard "lepton veto" of 45 GeV to reduce the background 
with very little effect on the signal. This cut slightly improves the reach of the analysis, but 
is most useful for increasing the purity (S/ B) of the Tevatron run I sample. 

The extra hard jets in the Wbb, Wcc, and single-top-quark backgrounds distinguish these 
processes from the signal. Since extra jets in the signal arise only in radiation from the hard 
process in 10-20% of the events, we consider a "jet veto" that removes any event containing 
a second jet with Exj > 30 GeV. The main effect of this cut is to improve the purity of the 
sample (S/B). The "jet veto" improves the significance (S/\fB) by a few percent at the 
Tevatron for high mx , but decreases it by a few percent at the LHC for all masses. This cut 
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is most useful if evidence of the signal is found and a higher purity confirmation is desired. 



B. Results 

To obtain numerical results, we adopt a minimal supergravity framework || to compute 
the masses of the top squark and its decay products, the mixing angle 6j, and the relevant 
branching fractions. Later, by varying the parameters of the model, we show that it is the 
top-squark mass itself that is the dominant variable in probing the A" couplings. 

We begin with common scalar and fermion masses of m = 100 GeV and mm = 150 
GeV, respectively, at the grand unified theory (GUT) scale. We choose a trilinear coupling 
Aq = —300 GeV and the ratio of the Higgs vacuum expectation values tan/5 = 4. The 
absolute value of the Higgs mass parameter fi is fixed by electroweak symmetry breaking 
and is assumed positive. Superpartner masses and decay widths are calculated with ISAJET 
7.50 [2^] and ISAWIG 1.1 After evolution from the GUT scale to the weak scale, for 



this set of parameters one obtains = 167 GeV, m^o = 53 GeV, m~± = 95 GeV, and 
sin^ = 0.8. In Fig. [7] we show the mass contours for the t\, the x5, and the Xi ■ A recent 
experimental lower bound of 80.5 GeV on the mass of the lightest chargino xt [23] is used to 
exclude regions (ex) of small m^j. Regions in which tachyonic particles would be generated, 
electroweak symmetry would not be broken, or in which the lightest supersymmetric particle 
(LSP) is not the lightest neutralino Xi are marked as excluded theoretically (th). 

The top-squark mass grows slowly with mo- The lightest chargino xt an d neutralino Xi 
masses are nearly independent of the common scalar mass uiq. In order to focus principally 
on top-squark mass dependence, we vary m and keep the other supersymmetric parameters 
fixed. Since the gaugino masses depend primarily on the choice of mi/2, variation of mo 
allows us to vary mj without an appreciable change in the masses of the decay products 
and only a slow rise of sin9j to 1. The results obtained in this section therefore depend 
primarily on the top-squark mass and on which decay channels are allowed kinematically. 

We reconstruct a longitudinally invariant mass M = M?+ fix-i where = \J m x + p^ x 
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is the transverse mass for the lepton— 6-jet system {m 2 x = PjjL, P£ = P{f + P/*), and fi T is 
the magnitude of the missing transverse energy for the event. The mass M has the useful 
feature that it produces a peak at vnq x — m^o + nib- Reconstruction of the peak provides 
a measure of the mass difference between the top squark and the lightest neutralino. We 
replace the mass definition of Ref. |JJ] because it is not longitudinally boost-invariant and 
hence is more sensitive to the detailed modeling of the decays. An added benefit is that the 
reach in mass with run I data is improved. 

In Fig. |^ we show a representative reconstructed mass M distribution at run II of the 
Tevatron (yS = 2 TeV) for a top squark of mass mt = 255 GeV. The coupling A" = 0.1 is 
chosen to be one order-of-magnitude below the current bound. The total background (B) is 
shown with its components: Wj includes Wc, Wj, Wbb, and Wcc; t includes all single-top- 
quark production modes. Both the "lepton veto" and "jet veto" are applied in this figure. 
The decrease in the background for M < 140 GeV evident in the figure is attributed to the 
cuts we impose. 

It should be easy to observe the sizeable deviation from the background distribution 
associated with a top squark produced with coupling strength A" = 0.1, shown in Fig. 
In Fig. [5] we show the distribution in M for the minimum value of A" = 0.04 required for 
a ha discovery with an integrated luminosity of 2 fb' 1 at \/~S = 2 TeV. The significance is 
calculated in a bin of width ±15% about the center of the peak in M. A bin width of ±10% or 
±25% produces the same significance to within a few percent. In this plot, — m^o ~ 201 
GeV. When this difference is reduced to ~ 150 GeV, the signal and background spectra 
peak in the same location, and sensitivity to the signal begins to be lost. 

Performing the same mass reconstruction for a series of top-squark masses, we examine 
quantitatively the functional dependence of Eq. (|5|) on A". For most of MSUGRA parameter 
space, aT~ p = 0.01, and bT^l = 0.5. The resulting significance (S/yfB) as a function of A" 
is shown in Figs. [10] and |11] for each run of the Tevatron. For A" < 0.2, the significance 
is proportional to {\") 2 /\fB. This behavior has two important consequences. In order to 
improve the bounds on the A" couplings by a further factor of 2 at run II, an additional 
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factor of 16 in integrated luminosity would be required. Thus, it will be difficult to set better 
limits than those presented here. On the other hand, the full run II integrated luminosity is 
not necessary to approach the limits. Second, the calculation is stable against uncertainties 
in the background estimation, or higher order corrections to the signal. Even a factor of 
2 uncertainty in the background would not shift the limits by more than 20%. If A" ~ 1, 
the loss of signal is important only for the highest masses we probe. This region might be 
covered by looking for a peak at large Exb in the many-jet plus 6-tag sample, where the 
chargino has decayed to three jets. 

In Fig. [12] we show the reach in A" for 165 < mj < 550 GeV for run I and run II of the 
Tevatron, and with the first 10 fb _1 at the LHC With an integrated luminosity of 2 fb _1 
at v^S* = 2 TeV, discovery at the level of 5a is possible provided that A" > 0.02-0.1 for 
top squarks with 165 < < 400 GeV. Otherwise, a 95% confidence-level exclusion limit 
can be set for A" as a function of the mass out to about 500 GeV. For the lower integrated 
luminosity and energy of the existing run I data, values of A" > 0.04-0.3 can be excluded 
at the 95% confidence level if mj = 165-350 GeV. We note that the run I results in Fig. |T2" 



represent a significant improvement over those in our Letter ||. The principal reason for 
this change is use of the new mass variable M. With the large increase in cross section at 
the LHC, both discovery regions and limits may be extended another order of magnitude in 
A", and to top-squark masses approaching 1 TeV. 

There are two reasons our analysis cannot be extended much below mj ~ 165 GeV. One 
is the loss of signal owing to cuts. While the cuts might be relaxed slightly, they are now 
close to the limit. The other is a loss of confidence in modeling of the background near and 
below the peak in the W — j mass distribution. Below this peak, mismeasurement problems 
and detector effects are significant sources of uncertainty. A detailed experimental analysis 
may be able to overcome this constraint. 

The virtue of probing the X'^ k couplings directly is that they may be bounded inde- 
pendently if no signal is found. The limits on individual couplings are determined by the 



relative weights of the corresponding parton luminosities. In Figs. |13| and [14], we show the 

14 



95% confidence-level limits that may be placed on A 312 , A 313 , and A 323 independently, as a 
function of top-squark mass at run I and II of the Tevatron. For = 255 GeV, absence 
of a signal in run I data means that A 312 < 0.08, and absence in run II data means that 
A 312 < 0.03. The bound on A 312 is strongest because its contribution is proportional to the 
d (g) s parton luminosity. Even at run I, a significant improvement can be made in the limits 
on A 312 and A 3 ' 13 . 

Another aspect of Eq. @ is apparent in the crossing of the limits for a common A" and 
A 312 . The three possible couplings contribute almost equally to the R p - violating decays of 
the top squark and chargino. In contrast, the production rate arises mostly from the A 312 
term. Hence, there is a net gain in measurable cross section in the i? p -conserving decay 
channel if only this coupling is large. At run II, a limit of A 312 < 0.2 can be placed for 
top-squark masses up to 550 GeV. Even A 323 may be probed directly if m? < 300 GeV. 

C. MSUGRA 

A specific set of parameters is used in the prior section to generate the masses of the 
top squark and its decay products. In this section we undertake a broader examination of 
the minimal supergravity model (MSUGRA) parameter space, and we present limits as a 
function of mo and m.1/2 for various choices of Aq, tan/3, and the sign of /i. This study 
allows us to conclude that for most of MSUGRA parameter space the top-squark discovery 
potential and bounds depend most strongly on the the top-squark mass alone. 

In Fig. |l^ we show contours in the mo — mi/2 plane of the minimum A" necessary for a 5a 
discovery of the top squark given 2 fb™ 1 of integrated luminosity at run II of the Tevatron. 
Plots are shown for Aq = —300 GeV, tan (3 = 4 and tan/3 = 30, and both signs of fi. The 
regions of experimental and theoretical exclusion are explained above, in conjunction with 
Fig. |7|. At low mo the rv\\ decay of xt saturates its i? p -conserving branching fraction. 
Naively, this effect would seem to make the electron and muon signal disappear. However, 
a small fraction of the taus satisfies the electron or muon tagging conditions. The net effect 
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is an increase in the mass reach in the region of small mo, most evident in the contour plots 
for tan /3 — 4. 



The discovery contours of Fig. 15 follow closely the top-squark mass contours. This 



point becomes evident if one compares Fig. [7] and the top-left-hand pane of Fig. [15] for 
which the same MSUGRA parameters are used. For example, the contour for A" = 0.1 in 
the top-left-pane of Fig. [15| can be seen to correspond to = 350-375 GeV in Fig. [7|. The 



same statement is true also for large tan/3 and both signs of \i: A" = 0.1 corresponds to 
m ~ = 350-375 GeV. 

The only region of significant dependence on the MSUGRA parameters, beyond those 
that set the top-squark mass, is at low tan (3 and /i < 0. There is a large suppression of the 
leptonic branching fractions of the chargino when |/i| ~ mi/2 ~ m . This suppression is due 
mostly to a destructive interference between the W boson and sneutrino (u) decay modes 



in this region po| . In terms of physical parameters, m^ is close to and the sign of the 
W — v interference term is the same as the sign of [i. Hence only the region of negative /i 
is affected. However, the cancellation decouples, as expected, as m^ becomes heavier than 
my/. Dependence of the bounds on the trilinear coupling A is negligible except in the region 
of small tan (3 and fi < where it serves as a moderator of the cancellation. As A becomes 
more positive, the cancellation increases. 

The results in Fig. ^ demonstrate that the reach in top-squark mass is significantly 
larger in single-top-squark production than it is for pair production, where the latter is 
restricted to the lower left-hand corner of the plot (m < 300 GeV, m^ < 170 GeV) ||24|| . 

Contours in the m — m.1/2 plane of the 95% confidence-level limits (1.96a) are shown in 
Fig. [16] for the same parameters as in Fig. [15|. Except for the same region of suppression at 
low tan f3 and fi < 0, limits can be improved with data from run II at the Tevatron by at 
least an order-of- magnitude for a broad range of m < 800 GeV and mi/ 2 < 300 GeV. The 
exclusion contours of Fig. [U] also follow closely the top-squark mass contours. The contour 
for A" = 0.1 in the top- left-pane of Fig. |H] can be seen to correspond to mr = 450 GeV 
in Fig. 0. This remains true for large tan/3 and both signs of [i: X" = 0.1 corresponds to 
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mr = 450 GeV. 

ti 

The strength of a direct probe of the _R p -violating couplings is that limits may be ex- 
tracted for each X'^- k coupling separately. In Fig. [T7|, we present contours in the m — r%/2 
plane of the 95% confidence-level limits that may be placed on A 3 ' 12 , A 3 ' 13 , or A323 on the 
assumption that only one coupling is non- negligible. The tan f3 = 4, /j, > limits for a 
common A" from Fig. [16] are also shown for comparison. Contours for each \'^- k follow con- 
tours of constant top-squark mass just as for a common A". Hence, the top-squark masses 
corresponding to A^, A^, or A323 = 0.1 are « 450 GeV, 350 GeV, or 215 GeV, re- 
spectively. In terms of a minimal supergravity model, the limits achievable at the Tevatron 
on individual \'L k couplings can be estimated for any value of m and mi/2 by comparing 
Fig. ^ and Fig. 0TJ or Fig. |T|. 



V. SUMMARY 

In this paper, we extend our study of s-channel production of a single top squark in 
hadron collisions through an R p - violating mechanism [|]]. Even if the R p - violating couplings 
A" are as small as 0.01, two orders of magnitude below the current bounds, the single-top- 
squark rate exceeds the pair rate for rrv^ > 100 GeV. This enhancement will allow searches 
at the Fermilab Tevatron to reach much higher top-squark masses than those based on 
traditional squark pair-production. 

Because of the overwhelming jet rate from standard strong interaction processes, R p - 
violating production of top squarks followed by _R p -violating decay into a pair of jets is not 
a viable means to exclude i?-parity-violating production of top squarks or to obtain useful 
bounds on mj or A". 

We study observability of single top squarks by focusing on the clean i? p -conserving 
decay, t\ —> 6%i , with xt l + + v i + Xi- The combined branching fraction for this cascade 
decay is of order 20% or greater for A" < 0.1, a region of significant interest. We address the 
case where the neutralino either decays outside of the detector volume or leaves remnants 
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that are too soft to identify; its signature is missing energy in the detector. We simulate 
both the signal and standard model background processes, including parton showering and 
hadronization as well as a full detector simulation. We investigate sensitivity to parameters 
in our analysis by exploring a wide range of the parameter space of the minimal supergravity 
model of supersymmetry breaking. 

With an integrated luminosity of 2 fb _1 at a/5* = 2 TeV, we show that discovery of a 
top squark with < 400 GeV at the level of 5a is possible provided that the i? p -violating 
coupling A" > 0.02-0.1. Otherwise, a 95% confidence-level exclusion limit can be set for A" as 
a function of the mass out to about 500 GeV. For the lower integrated luminosity and energy 
of the existing run I data, values of A" > 0.05-0.2 can be excluded at the 95% confidence 
level if = 165-350 GeV. These limits would constitute a significant improvement over 
the current 95% confidence-level upper bounds of A" < 1 for squarks of the third generation. 
With the large increase in cross section at the CERN Large Hadron Collider, the discovery 
region and limits may be improved by another order of magnitude in A" and extended 
to top-squark masses approaching 1 TeV. The discovery reach, or exclusion limits, at the 
Tevatron and LHC depend almost exclusively on the value of the top-squark mass, largely 
independent of other MSUGRA parameters. Since the i?-parity violation occurs only at the 
production stage, the X'^ k couplings may be bounded independently if no signal is found. 
Even at run I, a significant improvement can be made in the limits on Ag 12 and \ 3l3 . 

We propose a simple method to directly probe i?-parity-violating couplings at hadron 
colliders. This analysis searches for a single top squark that is produced in the s-channel 
via a baryon-number-violating interaction, and decays via standard i?-parity-conserving 
interactions. This strategy should be extended to other squarks and i?-parity-violating 
couplings. 
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TABLES 

TABLE I. Cuts used to simulate the acceptance of the detector for the Tevatron run II and 
run I (in parentheses if different). At the LHC, \rj\ < 2.5 for all objects. 

\Vb\ < 2 (1) E Tb > 40 GeV 

\r)i\ < 1.5 (1.1) E Tl > 15 GeV (20 GeV) 

M<2.5 ^Tj>30GeV 

E T > 20 GeV 
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FIGURES 




FIG. 1. Feynman diagram for the s-channel production of a single top squark via A3. 




m~ (GeV) 

FIG. 2. Cross section versus top-squark mass for i?-parity-violating production of a single 
top squark at run II of the Fermilab Tevatron (\^S = 2 TeV) with X'^ k = 0.1 compared with the 
i?-parity-conserving production cross section for top-squark pairs. 
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FIG. 3. Dijet mass distribution for the i?-parity-violating decay of the top squark at several 
masses (150, 200, 250, 300, 400, 500 GeV), compared with CDF dijet data from run I of the 
Tevatron jl~5| . 
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FIG. 4. Branching ratios for the top squark to decay into two jets and for the chargino xt to 
decay into three jets, both as a function of the i?-parity-violating coupling A". Also shown is the 
branching ratio for the top squark to decay into blvxi- 
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FIG. 5. The transverse energy Et spectrum (fb/GeV) of the tagged fr-jet for top squarks of 
mass 167 GeV and 402 GeV, and the Wbb background at the Tevatron \/~S = 2 TeV. We require 
Exb > 40 GeV, marked by the vertical line. 
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FIG. 6. The transverse momentum pr spectrum (fb/GeV) of the tagged lepton for top squarks 
of 167 GeV and 402 GeV, and the Wc background at the Tevatron = 2 TeV. We require 
Pti < 45 GeV marked by the vertical line. 



25 
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m (GeV) 

FIG. 7. Contours of constant top-squark mass m^, chargino mass m~+ and neutralino mass 



m^o (GeV) through the m - plane in MSUGRA with tan (3 = 4, A = -300 GeV, and ji > 0. 
Theoretically excluded regions are indicated with (th) and experimentally excluded regions with 
(ex). 
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FIG. 8. The reconstructed mass M distribution for single-top-squark production (S) and back- 
grounds (B) at the Tevatron (v 7 ^ = 2 TeV) for a top-squark mass = 255 GeV and coupling 
A" = 0.1 with all cuts applied. The Wj component of the background includes Wc, Wj, Wbb, and 
Wcc. The t component of the background includes all single-top-quark production modes (tbj is 
dominant; tb and Wt are negligible). 
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M (GeV) 

FIG. 9. The reconstructed mass M distribution for single-top-squark production (S) and back- 
grounds (B) at the Tevatron (y/S = 2 TeV) for a top-squark mass mj- = 255 GeV with all cuts 
applied. The coupling A" = 0.04 produces the minimum signal for a 5a significance at this mass. 
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FIG. 10. Statistical significance of the single-top-squark signal (S/y/B) in run I of the Tevatron 
(\J~S = 1.8 TeV, 110 pb _1 ) versus A" for a variety of top-squark masses (GeV). 
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FIG. 11. Statistical significance of the single-top-squark signal (S/ y/B) in run II of the Tevatron 
(VS = 2 TeV, 2 fb" 1 ) versus A" for a variety of top-squark masses (GeV). 
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FIG. 12. Lower limits on discovery (S/y/B = 5, solid), evidence (S/VB = 3, dashed), and 95% 
confidence- level exclusion {S/\fB = 1.96, dotted) for A" versus top-squark mass in run I of the 
Tevatron (y/S = 1.8 TeV, 110 ptr 1 ), run II of the Tevatron (y/S = 2 TeV, 2 fbr 1 ), and one year 
at the LHC (y/S = 14 TeV, 10 fb" 1 ). 
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FIG. 13. 95% confidence-level exclusion limits (S/y/B = 1.96) for \'^- k versus top-squark mass 
in run I of the Tevatron (y/S = 1.8 TeV, 110 pb" 1 ). The curve marked by A" is obtained if all 
three \'L k are set equal. There is insufficient rate to obtain a 95% limit on A323 at run I. 




FIG. 14. 95% confidence-level exclusion limits (S/y/B = 1.96) for Ag - fc versus top-squark mass 
in run II of the Tevatron (\/S = 2 TeV, 2 ftr 1 ). 
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FIG. 15. Contours of the minimum A" needed for a (5a) discovery with 2 fb _1 of data at run 
II of the Tevatron. Four mo — MSUGRA planes are shown with Aq = —300 GeV, \i > or 
[j, < 0, and tan/3 = 4 or tan/3 = 30. Theoretically excluded regions are indicated with (th), and 
experimentally excluded regions are indicated with (ex). For a given A", discovery should occur 
for values of tuq — m 1 / 2 below and left of the designated contour. 
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FIG. 16. 95% confidence-level (1.96<r) exclusion contours that can be set for A" with 2 fb _1 of 
data at run II of the Tevatron. Four mo — m 1 / 2 MSUGRA planes are shown with Aq = —300 GeV, 
[i > or [i < 0, and tan/3 = 4 or tan/3 = 30. Theoretically excluded regions are indicated with 
(th), and experimentally excluded regions are indicated with (ex). Absence of a signal excludes 
values of A" in the region below and left of the designated contour. 
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FIG. 17. 95% confidence-level (1.96<r) exclusion contours that can be set for X'^ 12 , or A323, 
with 2 fb _1 of data at run II of the Tevatron. Four too — m 1 /2 MSUGRA planes are shown with 
Aq = —300 GeV, n > 0, and tan/3 = 4. Theoretically excluded regions are indicated with (th), 
and experimentally excluded regions are indicated with (ex). Absence of a signal excludes values 
of the A3 jk couplings in the region below and left of the designated contour. 
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